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Abstract. - A novel collective excitation is predicted to exist at metal surfaces where a 
two-dimensional surface-state band coexists with the underlying three-dimensional continuum. 
This is a low-energy acoustic plasmon with linear dispersion at small wave vectors. Since new 
modern spectroscopies are especially sensitive to surface dynamics near the Fermi level, the 
existence of surface-state induced acoustic plasmons is expected to play a key role in a large 
variety of new phenomena and to create situations with potentially new physics. 



The long-range nature of the Coulomb interaction between electrons in a three-dimensional 
(3D) metal is known to yield collective behaviour, manifesting itself in the form of plasma 
oscillations [1]. These collective charge-density oscillations are a fundamental property of 
metals, and are known to be of basic importance in the interpretation of a large variety of 
experiments. 

In 1957 Ritchie introduced the idea of surface plasmons [2], which are collective oscillations 
bound to a metal surface. In the long-wavelength limit, the energy of the surface plasmon is 
given by huss = hijjp/\/2, where ujp = (47rne^/me)^''^ is the bulk-plasma frequency, n being 
the electron density and rUf. the electron mass. The concept of surface plasmons has played 
a significant role in a variety of areas of fundamental and applied research, from surface 
dynamics [3,4] to surface-plasmon microscopy [5], surface-plasmon resonance technology [6], 
and a wide range of photonic applications [7-9] . 

Another collective excitation arises in a two-dimensional (2D) gas of electrons confined to 
a plane. In this case, the 2D plasmon has energy h,uj2D = ?i(27rn^-^e^q/me)^^^ {v?^ is the 2D 
electron density), which goes to zero as q^/'^ as the 2D wave number q tends to zero [10]. The 
2D plasmon was first observed in artificially structured semiconductors [11] and more recently 
in a metallic surface-state band on a silicon surface [12]. 
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Since the typical energy huip of a three-dimensional (3D) plasmon is a few electronvolts, 
thermal excitation of bulk and surface plasmons is improbable and the electronic properties 
near the Fermi level cannot be influenced by these excitations. Plasmons in a 2D electron gas 
have lower energies, but they still do not affect electron-hole (e-h) and phonon dynamics near 
the Fermi level due to their square-root dependence on the wave vector. Much more effective 
than ordinary 3D or 2D plasmons in mediating, e.g., superconductivity would be the so-called 
acoustic plasmons with sound-like long-wavelength dispersion, which have spurred over the 
years a remarkable interest and research activity [13]. 

Pines originally suggested that acoustic plasmons could be realized in the collective motion 
of a system of two types of electronic carriers [14]. The possibility of having a longitudinal 
acoustic mode in a metal-insulator-semiconductor structure was anticipated by Chaplik [15]. 
Chaplik considered a simplified model in which a 2D electron gas is separated from a semi- 
infinite metal. He found that the screening of valence electrons in the metal changes the 
2D plasmon energy from the square-root behaviour to a linear dispersion, which was also 
discussed by Gumhalter [16] in his study of transient interactions of surface-state e-h pairs at 
metal surfaces. However, acoustic plasmons were only expected to exist for spatially separated 
plasmas, as pointed out by Das Sarma and Madhukar [17]. Acoustic plasma oscillations were 
observed in two-dimensionally confined and spatially separated multicomponent structures 
such as quantum wells and heterojunctions [18], and were then proposed as possible candidates 
to mediate the attractive interaction leading to the formation of Cooper pairs in high- Tc 
superconductors [19,20]. 

In this Letter, we show that in metals with a partially occupied surface-state band the 
dynamical screening at the surface provides a mechanism for the existence of a new acoustic 
collective mode whose energy exhibits a linear dependence on the 2D wave number q. We 
consider the example of the Be(OOOl) surface, with no d bands, which presents a very-highly 
populated two-dimensional Shockley surface-state band of the s-p^ symmetry [21-23], and 
as such is an ideal candidate to investigate the nature of the dynamical screening in a 2D 
electron gas that is immersed in a 3D system. The novel result of our detailed band-structure 
calculations, in which transitions between the surface-state band and all other states of the 
semi-infinite metal are fully included, is that in a real metal surface where a partially occupied 
quasi-2D surface-state band coexists in the same region of space with the underlying 3D 
continuum, a well-defined surface acoustic plasmon exists. 

First of all, we calculate the eigenfunctions and eigenvalues of a one-dimensional single- 
particle hamiltonian that describes the main features of the surface band structure [24,25]. 
Then we evaluate the dynamical density-response function x^{z,z';q,ui) of non-interacting 
electrons, and solve an integral equation to obtain the interacting density-response function 
X{z,z';q,u) in the random-phase approximation (RPA) which is known to be exact in the 
g' — > limit [26]. x{z, z'; q,u)) contains both bulk and surface states and accounts for all 
possible transitions between them. 

Assuming translational invariancc in the plane of the surface, the Fourier components of 
the induced density satisfy the equation 




(1) 



and the collective oscillations created by an external potential of the form 



r''\z;q,u) 



(27r/g)e' 
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Fig. 1 - RPA calculation of the imaginary part of the electron density induced in Be(OOOl), as 
obtained from eqs. (0 and for the wave number q = 0.05 versus the energy uj and the coordinate 
z normal to the surface. The crystal edge {z — 0) is chosen to be located half a lattice spacing 
beyond the last atomic layer, and 2: < corresponds to the interior of the solid, (a) 2D plot of 
lm[Sn{z; q, u!)]/uj . Positive and negative values are represented by red and blue colours, respectively, 
more intense colours corresponding to larger absolute values. The probability density of the Shockley 
surface state is represented on the left-hand side (red shaded areas), as a function of the 2 coordinate, 
(b) 3D plot of lm[Sn{z; q,Lj)]/uj. 



can be traced to the peaks of the imaginary part of the surface response function [27, 28] 

Ini[g(g,w)]= J dze''nm[Sn{z;q,Lj)]. (3) 

In fig. n we show our fuU RPA calculation of the imaginary part of the electron density 
Sn{z;q,u!) induced in a Be(OOOl) surface by the external electrostatic potential of eq. (0) 
with q = 0.05. This figure exhibits two distinct spectral features. The first feature occurs at 
the surface-plasmon frequency of valence (2s^) electrons in Be {hujs ~ 13 eV), where Im[5n] 
exhibits a strong maximum at the surface (z = 0) and pronounced Friedel oscillations in the 
interior of the solid (z < 0); these are the main characteristics of the ordinary monopolar 
surface-plasmon mode, which would also be present in a semi-infinite free-electron gas (FEG) 
of bulk density equal to that of valence electrons in Be. In addition, however, very localized 
local maxima are visible near the surface at an energy of 0.6 eV, which correspond to a low- 
energy collective excitation that would be absent in a semi-infinite FEG. A comparison with 
the probability density of the partially occupied Shockley surface state, which is shown by red 
shaded areas on the left-hand side of fig. QJi, clearly indicates that the low-energy collective 
excitation originates from this 2D surface-state band. Nevertheless, such a 2D electron gas 
alone would only support a plasmon that for q — 0.05 has energy uj2d = 2.7 eV, well above the 
low-energy excitation that is visible in fig. and it is only the combination of the strongly 
localized 2D surface-state band with 3D bulk states which allows the formation of this new 
mode. 

Figure |21 shows the imaginary part of the surface-response function g{q^Ld) of Be(OOOl), 
as obtained from eq. Q for increasing values of q. As follows from the figure, the excitation 
spectra is dominated by the conventional surface plasmon at hiOg ~ 13 eV, which can be traced 
to the characteristic pole that the surface response function g{q,uj) of a bounded 3D electron 
gas exhibits at this energy, and a well-defined low-energy peak with linear dispersion. Figure 
12 clearly shows that the low-frequency mode has linear dispersion, with the sound velocity 
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Fig. 2 - Energy-loss function lm[g{q,Lu)]/uj of Be(OOOl) versus the excitation energy u, as obtained 
from eq. Q for various values of the wave number: q = 0.05 (top panel), 0.1 (middle panel), and 0.15 
(bottom panel), in units of the inverse Bohr radius a^^. The peaks are dictated by the corresponding 
poles of the surface response function g[q,uj). In the long-wavelength limit (g 0), g{q,uj) is simply 
the total electron density induced by the potential of eq. 

being very close to the 2D Fermi velocity vjf. 

In fig. Owe show the energy of the low-energy mode in Be(OOOl) versus q (red solid line), as 
derived from the maxima of our calculated surface loss function lT[i[g{q,u;)], together with the 
well-defined plasmon energies that we obtain when only the surface-state band is considered in 
the evaluation of the dynamical density- response function (blue solid line). While the plasmon 
energies of electrons in the isolated surface-state band nicely reproduce in the long-wavelength 
limit the conventional plasmon dispersion lu2d of electrons in a 2D electron gas (green dashed 
line), the combination of this surface-state band with the underlying 3D system yields a new 
distinct mode whose energy lies just above the upper edge 

ujI^ ^ yjpq ^ ^2/2 of the 2D e-h 
pair continuum (shaded area). At small q, this surface acoustic mode has linear dispersion, 

7 -, 1 
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Fig. 3 - Dispersion of the new low-energy collective mode that is visible in fig. |21 (red solid line), 
together with the result we obtain when bulk states are omitted in the evaluation of the energy-loss 
function (blue solid line) and the conventional plasmon dispersion cjp^ = v'jF^ of electrons in a 
purely 2D electron gas with = 0.45 uq (green dashed line). The dotted line represents the upper 
edge tjj;^ = vWq -\- q^ /2 of the 2D e-h pair continuum (shaded area). We note that momentum and 
energy conservation prevents 2D e-h pairs to be produced for energies above lo^ ■ 
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with the sound velocity being very close to the Fermi velocity of the 2D surface-state band. 

Figures show that acoustic surface plasmons are low-energy collective oscillations 
appropriate to very localized slow Shockley electrons with Fermi velocity vjP < vq plus their 
associated cloud of swift bulk electrons with VpP ~ vq . These excitations cannot occur at the 
characteristic energy of the usual 2D plasmons, where Landau damping would be possible by 
virtue of 3D e-h pair creation. Instead, an acoustic plasmon mode occurs at energies that 
follow closely the upper edge of the 2D e-h pair continuum. 

In conclusion, we have shown that the combination of a localized surface-state 2D band 
of slow electrons with a bulk 3D electron gas of fast electrons {v^P > vjP), which takes place 
in the (0001) surface of Be, leads to the formation of a novel low-energy collective excitation. 
This new acoustic mode has linear dispersion at small wave vectors, with the sound velocity 
being very close to the Fermi velocity of the 2D surface-state band. A similar mode is also 
found to exist at the (111) surfaces of the noble metals [29]. We hope that our prediction of 
acoustic surface plasmons will stimulate new experiments in this direction. 
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